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' We develop a theory of collective phase description for globally coupled noisy excitable elements 

' exhibiting macroscopic oscillations. Collective phase equations describing macroscopic rhythms 

of the system are derived from Langevin-type equations of globally coupled active rotators via a 
O ' nonlinear Fokker-Planck equation. The theory is an extension of the conventional phase reduction 

' method for ordinary limit cycles to limit-cycle solutions in infinite-dimensional dynamical systems, 

such as the time-periodic solutions to nonlinear Fokker-Planck equations representing macroscopic 
rhythms. We demonstrate that the type of the collective phase sensitivity function near the onset 
of collective oscillations crucially depends on the type of the bifurcation, namely, it is type-1 for the 
saddle-node bifurcation and type-11 for the Hopf bifurcation. 
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I. INTRODUCTION 



^ ' The collective rhythms emerging from coupled dynamical elements are the most remarkable synchronization phe- 



nomena in nature [iHSl- These collective rhythms, or collective oscillations, are exhibited not only by populations 
of self-oscillatory elements, but also by populations of excitable elements, in which the elemental dynamics are not 
self-oscillatory but their coupling gives rise to nontrivial collective dynamics 'g'-^. The active rotator is a well-known, 
convenient phenomenological model that can describe both excitable and oscillatory dynamics. Systems of interacting 
active rotators have been extensively studied and their collective dynamics have been analyzed (ol-fTsj. 

Macroscopic synchronization between interacting groups of globally coupled phase oscillators exhibiting collective 
rhythms has attracted our attention, because many rhythms in the real world are actually collective rhythms generated 
by coupled networks of microscopic elements [l^, [l3| • The recent introduction of the Ott- Antonsen ansatz [l8l - [20| has 
facilitated theoretical investigations in this direction, although applicable only to a restricted class of coupled phase 
oscillators. As an alternative general approach to study the phase synchronization between macroscopic rhythms, we 
$-H have formulated a theory of collective phase description for coupled phase oscillators, which gives a reduced collective 
5^ 1 phase equation for the macroscopic rhythms [2l| - |2a | . With this method, we can analyze the synchronization properties 
between collective rhythms in a closed way at the macroscopic level. Each group of the oscillators can be treated 
as a single macroscopic oscillator, which facilitates theoretical investigations of interacting collective rhythms. In 
particular, we have derived a formula that relates the collective phase response of the macroscopic rhythm to the 
individual phase responses of the constituent oscillators [2l| - [25j . 

In this study, we develop a more generalized theory of collective phase description that is applicable to macroscopic 
rhythms generated by excitable elements as well as by oscillatory elements. Specifically, we consider interacting groups 
of globally coupled noisy active rotators and derive a coupled collective phase equation describing them. Our present 
theory can be applied to a broader class of interacting groups of globally coupled dynamical elements, and includes 
our previous theory for globally coupled noisy phase oscillators 21-23] as a special case. As an example, we analyze 
the type of the collective phase sensitivity function near the onset of collective oscillations, and reveal that it crucially 
depends on the type of bifurcation, as in the phase sensitivity functions of ordinary dynamical systems (26l. [27j. 
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Recently, Schwabedal and Pikovsky [28|, |29| studied somewhat similar issues on effective phase description of noise- 
induced oscillations in excitable systems. While their study treated a one-body problem for an ensemble of statistically 
independent excitable elements at the microscopic level, the present study focuses on globally coupled noisy excitable 
elements exhibiting collective oscillations, which can be considered a one-body problem at the macroscopic level. 
Specifically, the focus of this study is on the development of a phase reduction method for genuinely collective 
oscillations arising from mutual interactions in coupled dynamical systems. 

The organization of the present paper is as follows. In Sec. |TT1 we formulate a theory of collective phase description 
for globally coupled noisy excitable elements. In Sec. IIIIl we illustrate our theory using numerical simulations and 
reveal that the collective phase sensitivity function behaves very differently, depending on the bifurcation type of the 
collective oscillations. Concluding remarks are given in Sec. IIVI The attached appendices clarify the relation of our 
present theory to the conventional phase reduction theory of ordinary limit-cycle oscillators (App. present the 
results for a system of globally coupled phase oscillators with translational symmetry (App. [B]), and outline further 
extensions to general systems of globally coupled dynamical elements (App. [C|) . 

II. FORMULATION OF COLLECTIVE PHASE DESCRIPTION 

In this section, we formulate a theory of collective phase description for globally coupled noisy active rotators 
exhibiting macroscopic oscillations. The theory can be considered a phase reduction method for limit-cycle solutions 
in infinite-dimensional dynamical systems. The relation of the present theory to the conventional phase reduction 
of ordinary limit-cycle solutions is given in App. |^ The present theory is also an extension (which we explain in 
App. [B|) of our previous theory for phase oscillators plI - Esj . 

A. Langevin-type equations and the nonlinear Fokker-Planck equation 

We consider interacting groups of globally coupled active rotators described by the following equation: 

^ N N 

(t) = V {^f^) + 1 5: r 0^^) + v^ej^^ it) + e,z [^f^) (i) + 1 E ^i^O ' 

k=l k=l 

for j — 1, • • • ,N and (cr, r) = (1,2) or (2, 1), where 4>"^\t) is the state of the j-th active rotator at time t in the 
tJ-th group consisting of N elements. The first term on the right-hand side represents the intrinsic dynamics of the 
active rotator; the second term, the internal coupling between the elements in the same group; the third term, the 
noise independently applied to each rotator; the fourth term, the external forcing common to all rotators in the a-ih 
group; and the last term, the external coupling between the elements belonging to different groups. Each isolated 
active rotator obeys = «(</>), which describes either excitable or oscillatory dynamics [9]. The internal and external 

coupling functions are given by T {(j3^^\ (j)-^^ ) and T cT{<i}'j'\<i}'k'')^ respectively. The characteristic intensity of the 
internal coupling within a group is scaled to unity, whereas that of the external coupling between the groups is given 
by Cg > 0. The state-dependent sensitivity function of the active rotator is Z(0^'^^), which corresponds to the phase 
sensitivity function in the case of phase oscillators. We assume that the rotator variable (p takes values in [0, 27r] and 
that the above functions, i.e., w(0), Z{(f)), r(0, 0'), and ro.T-(0, 0'), are 27r-periodic in the respective state variables. 
The external forcing is denoted hy pa-{t), whose characteristic intensity is given by Cp > 0. The noise S,j'^\t) is assumed 
to be independent white Gaussian noise [s^, HH , whose statistics are given by 

(e]^' it)) = 0, (e]'^^ it) is)) = 2S,kS.rS it-s). (2) 

The noise intensity is characterized by _D > 0. 

In the continuum limit, i.e., N — cx), the Langevin-type equation ([T]) describing coupled groups of active rotators 
can be transformed into the coupled nonlinear Fokker-Planck equation as follows (see also Refs. [3, 0, IS, 16, 21, 23]). 
Let /'^"'^ {(/), t) represent the one-body probability density function of rotator state in the cr-th group, which is 
normalized as 



r"d0/(-) (0,0 = 1- 
Jo 



(3) 
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Because the mean- field theory exactly holds for the case of global coupling in the continuum limit 0] , we can average 
the internal and external coupling terms in Eq. ([T]) by the one-body probability density functions, /'•'^^ (</>, t) and 
f'^'^\(j),t), respectively. Then the Langevin-type equation ([T|) takes the form of a single- rotator equation as follows: 



(4) 



where 



<j 

+ CpZ (4^)) p. (t) + eg r dcf,' (^f \ 0') /^"^ (0', . 



(5) 



The single-rotator Langevin-type equation ^ can be transformed into the following coupled nonlinear Fokker-Planck 
equation: 



d_ 

di- 



(6) 



or, explicitly, 

d_ 

dt 



(^) 
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27r 



d^'T^r (0,0')/'"^ (0',O/^"^ (0,0 



(7) 



The first two terms on the right-hand side of the nonlinear Fokker-Planck equation ([T]) represent the internal dynamics 
of the cr-th group, the third term represents the external forcing applied to the cr-th group, and the last term represents 
the external coupling between the cr-th group and the r-th group. When the external forcing and the external coupling 
are absent, i.e., when Cp = = 0, each group of active rotators obeying Eq. ([7]) is assumed to exhibit stable collective 
oscillation, i.e., a stable time-periodic solution. We further assume that this situation persists even if Cp and/or Cg 
become slightly positive. 



B. Collectively oscillating solution and its Floquet eigenfunctions 

Let us assume ep — eg — and focus on a single group. The group index cr is dropped in this subsection. Then the 
nonlinear Fokker-Planck equation ([7]) with ep = = can be written in the following form: 



^0) 



2tt 



yr(0,0')/(0', 0/(0,0 



(8) 



We now assume the existence of a stable time-periodic solution to the nonlinear Fokker-Planck equation (|5]). It is 
known that such a solution exists in a certain model of globally coupled active rotators [9^, which we analyze in 
Sec. mil The stable time-periodic solution to the nonlinear Fokker-Planck equation can be described by 



/(0,t) = /o(0,e(t)), e{t)^n, 



(9) 



where Q and f2 are the collective phase and collective frequency, respectively. Inserting Eq. ([9]) into the nonlinear 
Fokker-Planck equation ([8]), we find that /o(0, 0) satisfies the following equation: 



2tt 



v{4>)+ d0'r(0,0')/o(0',e) /o(0,e) 



52 

^7J7?/o(0,0)- 



(10) 



Let u{(j),Q,t) represent a small disturbance to the unperturbed collectively oscillating solution, and we consider a 
slightly perturbed solution 



/ (0,t) = /o(0, e(t)) + u(0, 8(0,0 ■ 



(11) 
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Equation ([8]) is then linearized in m((/>, i.e. 



^u(0,e,o = i('/',e)u(0,e,t). 



Here, the hnear operator L(0, 9) is given by 

L(0,e)u(0,e) 

where 



J {((>,&) - n 



de 



(12) 



(13) 



2ir 



7/(0,9)/ d0'r(0,0')/o (</-', e) 
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/o(<^,9)/ d0'r(0,0')"(0',e) 



+ i?|^-(0,e). 



(14) 



Note that L{(j), 9) is time-periodic through 9, and therefore, Eq. is a Floquet-type system with a periodic hnear 
operator. Defining the inner product of 27r-periodic functions as 



1 



271- 



2tt 



u* ((/), 9) , u (</), 0)J = 2^ / I d(t>u* (0, 9) u ((/), 9) 



we introduce an adjoint operator of 1/(0, 9) by 

'u* (0, 9) , L (0, 9) u (0, 9)] = [l* (0, 9) 7.* (0, 9) , u (0, 9)' 

By partial integration, the adjoint operator L*{(f), 9) is explicitly given by 

d 



where 



L* {^,e)u* (0,9) 
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J* (0,9) + f7 
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u* (0,9), 



(15) 



(16) 



(17) 



J* (0, 9) u* (0, 9) = i; (0) —u* (0, 9) + 



27r 



d0'r(0,0')/o(0',9) 



ru* (0,9) 



+ / d0' r (0', 0) /o (0', 9) — u* (0', 9) + D—u* (0, 9) 



(18) 



In the calculation below, we use the Floquet eigenfunctions of the periodic linear system Eq. (|12p and its adjoint 
system associated with the zero eigenvalue, i.e.. 



L(0,9)uo ('/',e) 
L* {^,e)u* (0,9) = 



J* (0,9) + f) 
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uo (0,9) =0, 
u* (0,9) =0. 



(19) 
(20) 

(21) 

as is confirmed by differentiating Eq. (fTO]) with respect to 9. Using the inner product (fTS]) with the right zero 
eigenfunction (pTjl . the left zero eigenfunction Uo(0, 9) is normalized as 



Note that the right zero eigenfunction mo(0j ©) can be chosen as 

"0 = ^/o('^,0) 



u* (0,9),uo (0,e) 



1 
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)wS(0,9)7/o(0,e) = i. 



(22) 
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Here, we should note that the foUowing equation holds (see also Eq. (jA16p and Ref. [3; 



d_ 



2tt 



d4>ui ((/.,e)uo(0,e) 



2n 



1 

0. 



d0 

2tt 



< (<^' 0) ae"" ('^' 0) + "0 ('^' 0) QQ< ('^' 0) 



(<^, e) J (0, 6) uo (0, e) - (0, e) r {4>, e) (0, e) 



Therefore, the following normalization condition is satisfied independently for each Q as 

dcf,u* (0,6) Mo (0,6)^1. 



(23) 



(24) 



In the following subsections, using the collectively oscillating solution and its Floquet zero eigenfunctions, we formulate 
a theory of collective phase description for globally coupled active rotators. 



C. Collective phase reduction 



Let us introduce the effects of weak external forcing and external coupling, i.e., we assume < ^ 1 and 
< Eg ^ 1. We treat the third term and the last term in Eq. ([7]) as perturbations. Using the idea of phase 
reduction _25, we can derive a coupled collective phase equation from the nonlinear Fokker-Planck equation ([T]). 
Namely, we project the dynamics of the nonlinear Fokker-Planck equation ([7]) onto the unperturbed collectively 
oscillating solution as 



d^u* (0,e(^)) — [z(0)/o(0,e("' 



P. (t) 



dci^u* (0, e('^)) ^ J^^ d4>' r,, (0, 4>') fo (0', e(^)) /o (0, e^^)) 



where we approximate f'^"\(j),t) by the unperturbed solution /o(0, 6^"^^) and use the fact that 

d 



d07.5(0,e)-/o(0,e) d^u*i^,e)uo{^,e) = n. 

Thus, the collective phase equation is obtained in the following form: 

e(^) {t) = n + (e(^)) {t) + e.r.. (e(^), e^^)) . 

Here, the collective phase sensitivity function is given by 



and the effective phase coupling function is given by 



271 



2-77 



Far (e("\e(^)) = deb dcb'r,r{<l>,ci)')ko (0,e(")) /o (0',e(^)) , 

where the kernel function is defined as 

fco (0,e)-/o(0,e)^u5(0,e). 

In the next subsection, we simplify the collective phase equation by applying the near-identity transformation. 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



6 



D. Normal form of the collective phase equation 

Applying the near- identity transformation 0, [H, [s^l to Eq. p7|. namely, by slightly (i.e., on the order of 0(eg)) 
changing the definition of the collective phase Q, we can obtain the following collective phase equation 

e('^)(t) = n + epC (e^'^)) p. {t) + eg^^r (e^^^ - e^^') , (31) 

where the collective phase coupling function is given by 




or, more explicitly, by 

I /'27r /'27r p'Itt 

7ar(e) = — / dx d4> d(/)'r,,(<^,(/)')fto(0,A + e)/o(0',A). (33) 

Jo Jo Jo 

Note that the derivation of Eq. pip from Eq. (I?fl) follows procedures developed in the appendix of Ref. [1^ where the 
details of the near-identity transformation are given. In addition, when the external forcing is absent, i.e., when = 0, 
Eq. ([5^ can be understood as the result of the averaging method [1] . In summary, the collective phase equation ([31]) 
with the collective phase sensitivity and the collective phase coupling functions, i.e., Eq. ([28)) and Eq. ([33)) . are derived 
from the Langevin-type equation ([l} via the nonlinear Fokker-Planck equation ([7]) . 



E. Adjoint method for the left zero eigenfunction 

Finally, we explain a numerical method for obtaining the left zero eigenfunction. From Eq. (|20p . the left zero 
eigenfunction Mg((/), 8) satisfies 

= (34) 

which can be transformed into 

d 

—ul (0, -Vts) = ,/* (0, -17s) ul ((/), -VLs) (35) 
OS 

with = — f2s. A relaxation method using Eq. (I55)) . which numerically calculates the eigenfunction associated with 
the zero eigenvalue by evolving Eq. (1351) until all eigenfunctions with negative eigenvalues decay, is a convenient 
method for obtaining the left zero eigenfunction (see also Eq. (IA23P and Refs. [27l. |32|. Issj ) . 

As mentioned at the beginning of this section, in App.lXJ we clarify the relation between our present formulation and 
the conventional phase reduction of ordinary limit-cycle solutions. In particular, it should be noted that the variable 
of the nonlinear Fokker-Planck equation plays the role of the vector component index in App. |X| In App. [B] it is 
shown that the present formulation includes our previous theory for phase oscillators [2l| - [23| as a special case with 
continuous translational symmetry in 0. 



III. NUMERICAL ANALYSIS OF GLOBALLY COUPLED ACTIVE ROTATORS 

In this section, we illustrate our theory (developed in the preceding section) using numerical simulations. First, 
we introduce a complex order parameter and its related quantities, which we will use to characterize the collective 
oscillations. Thereafter, we present the results of numerical simulations. In particular, we analyze the type of the 
collective phase sensitivity function near the onset of collective oscillations. 

A. Order parameters 

To quantify the collective oscillations, we introduce a complex order parameter 

R^t)e^m ^ Y.e^'^-'-'^ ^ d0e^^/(0,t), (36) 
k=i -^0 



7 



which is time-dependent. We also use the following real constant order parameter [P] 



(37) 



where (•)( is a long-time average. The order parameter S is non-zero only when collective oscillations exist. In 
addition, using the real and imaginary parts of the complex order parameter, we introduce the following pair of real 
order parameters: 



X = _Rcos<l>, 



Y = i?sin$. 



(38) 



In visualizing the limit-cycle orbit of the nonlinear Fokker-Planck equation with the infinite-dimensional state space, 
we project the collectively oscillating solution onto the X-Y plane as 



Xo(e) + iro(e) = i?o(e)e**«(®) = / dcbe'^fo{ci>,Q). 

JQ 



(39) 



B. Collectively oscillating solution and other associated functions 

In numerical simulations, we use the following model of globally coupled active rotators Q: 

V {(/))— u! — r siii{4)), w,r>0, 

and 



(40) 



r (^, <^') - r,, (<^, 0') = - sin (0 - 0') . (41) 

For simplicity, the state-dependent sensitivity function of the active rotator is assumed to be constant, i.e., 

Z(0) = 1. (42) 

Here, we summarize several properties of the above model Q : (i) The active rotators are excitable when uj < r, whereas 
they are oscillatory when uj > r. (ii) In the noiseless case, i.e., D = 0, the model exhibits collective oscillations under 
the condition w > r, which emerges via a saddle-node bifurcation, (iii) In the phase-oscillator limit, i.e., r = 0, 
the model exhibits collective oscillations under the condition D < Dc — 1/2, which occurs via a supercritical Hopf 
bifurcation. 

In our numerical simulations, the parameters of the active rotators are fixed at w = 0.50 and r = 0.52, values 
that satisfy the excitable condition lu < r. The noise intensity is mostly chosen to be D = 0.22, for which collective 
oscillations exist with the collective frequency ~ 0.1955. In numerical simulations of the nonlinear Fokker-Planck 
equation ([7]) and the adjoint equation ([35l) . we used the pseudospectral method with 2^" modes. 

Figure [1] shows the limit-cycle orbit projected on the X-Y plane, which was obtained from the nonlinear Fokker- 
Planck equation ([7]) with ep = Cg = 0, i.e., Eq. ([5]). The collectively oscillating solution /o(0, 6) and other associated 
functions, i.e., uo((/), 9), Uq((/), 9), fco(<^, O), and r(6,6'), are displayed in Fig. [2] 

The collective phase sensitivity function C(0) the collective phase coupling function "/^^[Q) are shown in 
Fig.[3Ka) and Fig.[3Ub), respectively. The collective phase sensitivity function ({Q) is quite different from the individual 
state-dependent sensitivity function Z{(f)) = 1. Despite the fact that Z{(j)) is constant, the collective phase sensitivity 
function depends on 8, in sharp contrast to the case of globally coupled phase oscillators where C(0) never 

depends on 6 if Z{(j)) is constant 22] (see also App. [B|). Consequently, although the external coupling function has 
only a first harmonic, i.e., rcrr(0, 0') = — sin(0 — <j)'), the collective phase coupling function 7crr(9) includes higher 
harmonics (and also a uniform component); that is also in contrast with the case of globally coupled phase oscillators 
in which jcrriQ) is sinusoidal whenever rcrT{4>) is sinusoidal ^23,] (see also App. [B|) . 

The collective phase sensitivity function C(0)) measured by direct numerical simulations of the system by applying 
sufficiently weak stimulating impulses with = 0.01, is compared to the theoretical curve in Fig. [31[c). The results of 
the nonlinear Fokker-Planck equation corresponding to the continuum limit N ^ oo agree perfectly with our theory; 
the results of the Langevin-type simulation with N = 10^ active rotators averaged over 100 ensembles agree with 
our theory within the fiuctuations due to finite size effects. Therefore, the formulation developed in Sec. Ullhas been 
validated by the above-described numerical simulations. 
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C. Types of collective phase sensitivity function 



Now, we investigate the dependence of the cohective phase sensitivity function C(0) on the noise intensity D. In 
particular, we analyze the type of the collective phase sensitivity function C(0) near the onset of collective oscillations. 

FigurelUa) shows the dependence of both the order parameter S and the collective frequency on the noise intensity 
D (see also Ref. Q). As the the noise intensity is increased, collective oscillations arise via a saddle- node bifurcation 
and disappear via a supercritical Hopf bifurcation, in consistent with the results found in Ref. In the present 
study, we further investigate the dependence of the collective phase sensitivity function on the noise intensity. The 
typical shapes of the collective phase sensitivity function ^(0) near the onset of collective oscillations are displayed 
in Fig. m^b). The collective phase sensitivity function near the saddle-node bifurcation {D — 0.10) satisfies C(0) > 0, 
which is referred to as type-I (26 . i27i] . That is, the collective phase is always advanced by positive weak impulses. 
In contrast, the collective phase sensitivity function near the Hopf bifurcation [D = 0.25) is type-II [IHUS]) that is, 
C(0) has negative parts as well as positive parts. 

The above results are the same as those for the finite-dimensional dynamical system [27j . However, note that 
they are the macroscopic collective characteristics of an infinite-dimensional dynamical system that emerged from 
interactions of individual, microscopic elements. 

As shown in this example of globally coupled active rotators, a coupled dynamical system may have multiple routes 
(such as saddle-node bifurcation or Hopf bifurcation) to collective oscillations depending on the system parameters, 
yielding different collective phase sensitivity functions. It is known that in weakly interacting ordinary limit-cycle 
oscillators, the difference in the type of phase sensitivity function can lead to very different synchronization dynam- 
ics ^2^ [23 • Therefore, the present results imply that when we consider populations of coupled dynamical elements 
undergoing collective oscillations, the bifurcation of each population leading to the collective oscillations can strongly 
influence the synchronization properties. 



We developed a theory of collective phase description for globally coupled noisy active rotators where macroscopic 
rhythms are generated by excitable elements. By projecting the dynamics of the nonlinear Fokker-Planck equation 
onto its stable limit-cycle orbit in the infinite-dimensional state space, we derived a collective phase equation describing 
macroscopic rhythms of the system. On the basis of our theory, we analyzed the type of the collective phase sensitivity 
function near the onset of collective oscillations. We obtained two distinct types of the collective phase sensitivity 
function; they correspond to two different bifurcation routes to collective oscillations. 

The theory we developed in this study can be considered a phase reduction method for limit-cycle solutions to 
infinite-dimensional dynamical systems (see App. and includes our previous theory for phase oscillators as a 
special case (see App.[B]). Furthermore, the extension of this approach to general systems of globally coupled dynamical 
elements is briefiy discussed in App. [Cj 

Finally, we remark the relation between the phase and the breaking of continuously translational symmetry, which is 
the core of the phase reduction theory . The phase of the ordinary limit-cycle solution is associated with temporal 



translational symmetry breaking. In Refs. j2l| - |23| . the collective phase of the nonlinear Fokker-Planck equation 



of globally coupled phase oscillators is assigned to a mixture of the spatial and temporal translational symmetry 



breakings, where the term "spatial" refers to the variable 4>. The formulation in Refs. |2l| - |23l | is also essentially the same 



as the phase dynamics of wavefronts in reaction-diffusion systems [2]. However, as in the case of ordinary limit cycles, 
spatial translational symmetry is not essential (although it is typically assumed in the wavefront dynamics of reaction- 
diffusion systems) for the phase description of limit-cycle solutions in infinite-dimensional dynamical systems. In this 
paper, the collective phase of the nonlinear Fokker-Planck equation is assigned only to the temporal translational 
symmetry breaking. Therefore, the formulation is applicable to nonlinear Fokker-Planck equations describing not 
only phase oscillators but also active rotators, even though the latter case does not possess spatial translational 
symmetry. A similar formulation for stable time-periodic solutions to reaction-diffusion systems will be presented in 



IV. CONCLUDING REMARKS 
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Appendix A: Conventional phase reduction of limit-cycle oscillators 

We briefly review the conventional phase reduction method for coupled limit-cycle oscillators @, [3^ [3^ to make a 
comparison to our present theory. This will clarify that the formulation in Sec.|TT]can be considered a phase reduction 
method for limit-cycle solutions in infinite-dimensional dynamical systems. 



1. Coupled limit-cycle oscillators 



We consider coupled limit-cycle oscillators described by the following model: 



X'"^ {t)=F[X 



(Al) 



for (ct, r) = (1, 2) or (2, 1), where X^*^-* (t) is the d-dimensional state of the cr-th limit-cycle oscillator at time t. The first 
term on the right-hand side represents the intrinsic dynamics of the limit-cycle oscillator, the second term represents 
the external forcing, and the last term represents the mutual coupling between the oscillators. The dynamics of the 
isolated limit-cycle oscillator are given by X = F{X). The external forcing is denoted by p„{t), whose characteristic 
intensity is given by tp > 0. The mutual coupling function is Go-r (^''^' , ^'■^■'), whose characteristic intensity is 
e,>0. 

2. Limit-cycle solution and its Floquet eigenvectors 

Let us assume tp — tg — and focus on a single isolated oscillator. The oscillator index a is dropped in this 
subsection. The limit-cycle solution to Eq. (jAll) without external forcing or mutual coupling can be described by 



(A2) 



where 9 and uj are phase and frequency, respectively. Inserting Eq. (jA2p into Eq. (|A1[) with ep = eg = 0, we find that 
Xf){6) satisfies the following equation: 



uj-Xo {9) = F{Xo {9)). 



(A3) 



Let u{9,t) represent a small disturbance to the limit-cycle solution; let us then consider a slightly perturbed solution 

X{t)^Xo{9)+ui9,t). (A4) 
Equation (jAl[) with = = is then linearized in u(9,t) as 



^^u{9,t)^L{9)u{9,t), 



where the linear operator L{9) is given by 



L (9) u (9) = 



with a Jacobi matrix 



J{9) = 



dF{Xoi9)) 
dXo{9) ■ 



ui9) 



Defining the inner product as 

\u* (9) , u {9) 
we introduce an adjoint operator of L{9) by 



2^ 7o 



d9u* [9) ■u[9), 



(A5) 



(A6) 



(A7) 



(A8) 



u* {9), L{9)u{9) = L* {9)u* {9),u{9) 



(A9) 
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The adjoint operator L*{9) is explicitly given by 

L* {e)u* [9) 

with the transposed matrix of the Jacobi matrix 



u* (9) 



J* (9) = J{9f 



(AlO) 



(All) 



We use the Floquet eigenvectors of the periodic linear system Eq. (|A5[) and its adjoint system associated with the 
zero eigenvalue, i.e., 



L{9) uo{9) 
L* {9)u* {9)^ 



Ji9)-LU 

J* {9)+uj 



uo (9) = 0, 
u* (9) = 0. 



Note that the right zero eigenvector Uo{9) can be chosen as 



Uo {9) = -Xo (9) , 



(A12) 
(A13) 

(A14) 



which is confirmed by differentiating Eq. (|A3|) with respect to 9. Using the inner product (|A8[) with the right zero 
eigenvector (|A14|) . the left zero eigenvector Uq{9) is normalized as 



«o (0) , Uo (9) 



2^ Jo 



d9u*o (9) ■uo{9) = 1. 



(A15) 



As in the infinite-dimensional case treated in the main text, the following equation holds [32 



d_ 
d9 



u*i9).uoi9) =«S(^)-^ + ^.noW 



1 r 
0. 



d9 d9 
u* {9) ■ J (9) uo {9) ~ J* (9) u*o (9) ■ uo {9) 



Therefore, the normalization condition is satisfied separately for each 9 as follows 

u* (9) ■ Uo {9) = 1. 

The left zero eigenvector is the so-called phase sensitivity function 0, i.e., 

u*o {9) = Z{9). 



(A16) 
(A17) 
(A18) 



3. Phase reduction and the normal form 

Now, we introduce the effects of weak external forcing and mutual coupling, and treat the second and last terms in 
Eq. (|Aip as perturbations. Using the phase reduction method [2], we can derive a phase equation from Eq. (jAl|) as 
follows: 



^('^) {t)^uj + epZ (0('^)J • {t) + tgZ j • G^r [Xo , Xo [9^^^)) . (A19) 

Applying the near- identity transformation 0, [25l [3^ , we can obtain the following simplified phase equation: 

^)(-) {t)=u + epZ (^?(-)) • p„ {t) + egV^r (e^^^ - 0^^^) , (A20) 
where the phase coupling function is given by 



r.r W - 7^ / dXZ{X + 9)- G„r {Xo{\ + 9),Xo{X)) 



(A21) 
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4. Adjoint method for the left zero eigenvector 

From Eq. (|A13|) . the left zero eigenvector Uq{9) satisfies 

u;-^u*{0)=-.ri0)u*ie), (A22) 

which can be transformed into 

-^u*n (~ujs) = J* (-ujs) u*n (-ujs) (A23) 
as 

with 9 = —UJS. As is well known [23, IH, [s^, the relaxation method using Eq. (IA23|) is a convenient method for 
obtaining the phase sensitivity ftmction Z(9) numerically. 



5. Remarks on the relation to the conventional phase reduction of ordinary limit cycles 

Now, we see that the formulation in Sec. |lT] clearly corresponds to that in this appendix. In particular, the following 
relations are important: 

fo{cf>,Q)^Xo{9), (A24) 

2-77 

d^u* (0, 6) uo (0, 8) ^ u*{9) ■ uo{9). (A25) 

Namely, the formulation in Sec.|ll]is considered a phase reduction of the limit-cycle solution to the infinite-dimensional 
dynamical system that is continuously parameterized by the variable (p. 

Appendix B: Phase oscillators with translational symmetry 

Here we clarify the relation between the formulation developed in Sec.|lT]and our previous formulation for globally 
coupled phase oscillators (2l| - |23| . Phase oscillators can be described as a special case of active rotators as follows: 

V (0,) -> (Bl) 

r(0„0fe)^r(0j-(/.fe). (B2) 

In this case, the nonlinear Fokker-Planck equation has a translational symmetry with respect to 4>. Reflecting this 
symmetry, the collectively oscillating solution and other associated functions have the following properties: 



/o( 




= /o(</ 


;.-9,o) 


= /o ['^) , 


Wo ( 


0,9) 


= Uo (' 


^-9,0) 


= Uo (95) = 


Uo ( 


4>,Q) 


= K) (' 


?!>- 9,0) 


= K) if ) , 


ko{ 


0,9) 


= fco (<; 


6-9,0) 


= fco (if) = 



(B3) 
(B4) 
(B5) 
(B6) 



^fo if) , 
dip 

• if) if) . 

In addition, the kernel function k{){ip) is normalized as 

difka{ip)^l d(p fo {(fi) —Uq (ip) = dipu^ (ip) Uq (ip) ^ 1. (B7) 
Jo dip Jo 

The collective phase sensitivity function C(9) is given by 

C(9)= / d0Z(0)A:o(0-9) = / dpZ{p + e)koip), (B8) 
Jo Jo 

whose 9-dependence exists only when the phase sensitivity function Z[(l>) is not constant. Furthermore, the external 
coupling function Fq-t of the phase oscillators can be written in the following form: 

r<.r(0f\4^^)^r..(0f'-0i^^). (B9) 
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Again, reflecting the translational symmetry, the effective phase coupling function F ar depends only on the collective 
phase difference and coincides with the collective phase coupling function 7o.r as follows: 



Fc 



(^0(-)^ 0(r)^ = F^r (e(") - e(^), O) = Far (Q^"^ - e'"') = 7ar (O^^^ - e(")) , (BIO) 



or, explicitly. 



2ir /.27r 



V Tar + e('^) - e^^) j fco iv) fo iv') . (Bii) 



Here, we should note that when the external coupling function T^t is sinusoidal, the collective phase coupling function 
•Jar also has a sinusoidal form. The above collective phase sensitivity and coupling functions for the phase oscillators, 
i.e., Eq. (|B8p and Eq. (|Blip . are equivalent to the formulae that we derived in Refs. [13, H^- Precisely speaking, the 
above system of notation is slightly different from the previous one [2l| - [23j as follows: 



fo M = fo i^) , (B12) 

uo if) = -uq (lp) , (B13) 

< iv) = -< iv>) , (B14) 

fco iv) = fo (^) (^) = -fo {^) -^K [ip) = fc* (^) . (B15) 



Here, the previous notations of Refs. |2l| - |23l | are denoted by foiv), uoi'p), ""o ^^"^ fco (</')• These differences, which 
are not essential, come from the fact that the right zero eigenfunction was chosen as uoi^) = dfo{^)/d(p, in contrast 
to Eq. (HH), i.e., uo{ip) = ~dfo{(p)/dip. 

Let us consider the following simple model: 

r(0) =r,,(,^) = -sin(0 + a), |a|<|- (B16) 

In this case, collective oscillations exist under the condition D < Dc ^ cos(q;)/2 0, H^. The collectively oscillating 
solution /o((/), 8) and other associated functions, i.e., mo(0, 8), UQ{(j>,Q), fco((?I>, 8), and F(8,8'), are shown in Fig. [5l 
where the parameters are set to a = Stt/S and D = Dc/2 — cos(q!)/4. All functions are translationally symmetric in 
sharp contrast to Fig. [5] Using the translational symmetry. Figs. [3a)(b)(c)(d) can be compressed into Fig. [HI which 
is essentially the same as the results in Ref. (23| . 



Appendix C: Extension to general systems of globally coupled elements 



The formulation in Sec.|TT]can be straightforwardly extended to general systems of globally coupled noisy elements, 
whose outline is sketched below. We can consider general dynamical elements such as excitable, oscillatory, or chaotic 
units, as long as the system exhibits collective oscillations described by a time-periodic solution. We consider a group 
of globally coupled general dynamical elements subjected to independent white Gaussian noise described by 



1 ^ 



k=l 



which can be transformed into the following nonlinear Fokkcr-Planck equation: 



d_ 

di 



f{X,t) 



d 
dX 



F{X)+J dX' G (X, X') / {X\ t) I / (X, t) 



D 



Therefore, a collectively oscillating solution can be described by 

= /o(x,e), Q{t) = n, 

whose right zero eigenfunction can be chosen as 

(X,8) = ^/o(X,8), 



d d 
dX ' dX 



f{X,t). 



(CI) 



(C2) 



(C3) 



(C4) 
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and whose left zero eigenfunction is normalized as 

,{x,e)uo{x,e) = i. (C5) 



j dXu*\ 



It is clear that collective phase description for general systems of globally coupled noisy elements exhibiting macro- 
scopic rhythms can be formulated, in principle, in the same way. 
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FIG. 1: (Color online) Globally coupled active rotators (AR) in the excitable regime exhibiting macroscopic rhythms, (a) 
Macroscopic limit-cycle orbit projected on the X-Y plane, (b) Wave forms of Xo(0) and Yo[Q). Parameters are uj = 0.50, 
r = 0.52, and D = 0.22, which give Q. ~ 0.1955. 
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FIG. 2: (Color online) Globally coupled active rotators (AR). (a) Collectively oscillating solution fo{(j),0). (b) Right zero 
eigenfunction uo((;/>,0). (c) Left zero eigenfunction uq{<I),Q). (d) Kernel function fco(<?!', O). (e) Effective phase coupling 
function F{Q, O ). Parameters are oj = 0.50, r = 0.52, and D = 0.22. 
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FIG. 3; (Color online) Globally coupled active rotators (AR). (a) Collective phase sensitivity function C(©)- (b) Collective 
phase coupling function 7(6). (c) Comparison of collective phase sensitivity functions C(Q) between theoretical curve (Theory) 
and direct numerical simulations with impulse intensity = 0.01, i.e., Langevin-type equation (Langevin) and nonlinear 
Fokker-Planck equation (NLFPE). Parameters are it: = 0.50, r = 0.52, and D = 0.22. In Langevin-type simulations, the 
number of active rotators in the system is A'^ = 10^ and the results are averaged over 100 ensembles, whose means and standard 
deviations are shown. 




FIG. 4: (Color online) Globally coupled active rotators (AR). (a) Dependence of order parameter S and collective frequency 
f2 on noise intensity D. (b) Collective phase sensitivity function C(©)- Parameters are l>j — 0.50 and r — 0.52. The collective 
phase sensitivity function near the saddle-node bifurcation is type-I, whereas that near the Hopf bifurcation is type-II. 
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FIG. 5: (Color online) Globally coupled phase oscillators (PO) corresponding to the translationally symmetric case, (a) 
Collectively oscillating solution /o(0, 0). (b) Right zero eigenfunction uo{(j),Q). (c) Left zero eigenfunction uq{(I),Q). (d) 
Kernel function ko{4i, O). (e) Effective phase coupling function F{0, O')- Parameters are a = Stt/S and D = Dc/2 = cos(a)/4. 
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FIG. 6: (Color online) Globally coupled phase oscillators (PO). (a) Collectively oscillating solution foif)- (b) Right zero 
eigenfunction uo{lp). (c) Left zero eigenfunction uq{lp). (d) Kernel function ko{ip). Parameters are a — 3n/8 and D — Dc/2 — 
cos(a)/4. 



